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Abstract 

CIZA J1358.9-4750 is a nearby {z = 0.074) pair of clusters of galaxies located close 
to the Galactic plane. It consists of two X-Ray extended humps at north-west and south¬ 
east separated by 14 arcmin (~ 1.2 Mpc), and an X-Ray bright bridge-like structure in 
between. With Suzaku, the south-east hump was shown to have a temperature of 5.6 ± 

0.2 keV and the north-west one 4.6 ± 0.2 keV. Neither humps exhibit significant central 
cool component. The bridge region has a temperature higher than 9 keV at the maximum, 
and this hot region is distributed almost orthogonal to the bridge axis in agreement with the 
shock heating seen in numerical simulations at an early phase of a head-on major merger. 

This resemblance is supported by good positional coincidence between the X-Ray peaks 
and cD galaxies associated with each cluster. In a short exposure XMM-Newton image, a 
significant intensity jump was found at a position where the Suzaku-measured temperature 
exhibits a steep gradient. These properties indicate the presence of a shock discontinuity. 

The Mach number is estimated to be 1.32 ± 0.22 from the temperature difference across 
the identified shock front, which gives the colliding velocity of approximately 1800 km 
s“^. From optical redshifts of the member galaxies, the two clusters are indicated to be 
merging nearly on the sky plane. Thus, CIZA 11358.9-4750 is considered as a valuable 
nearby example of early-phase merger with a clear shock feature. 

Key words: clusters of galaxiesi — merging clusters 2 — shock waves 
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1. Introdcution 


Galaxy clusters are eonsidered to grow through their mergers. Huge gravitational energy up 
to ~ 10®"^ erg is released during the event, espeeially in ease of “major merger”, i.e, nearly head- 
on merging of two elusters with similar sizes (e.g. Rieker and Sarazin 2001). Part of this energy is 
converted into thermal energy of the intraeluster medium (ICM), mainly via eollision-indueed shoeks. 
These shoeks also convert the energy into non-thermal energy sueh as those to aeeelerate eosmic-rays, 
drive ICM turbulence, and amplify intraeluster magnetie fields. However, it is still unknown how the 
released dynamieal energy is distributed into these different ehannels. 

One of the approaehes to tackle this problem is to observe merging elusters with elearly iden¬ 
tified eollision parameters, sueh as eolliding veloeity, viewing angle, and the original mass of the 
eolliding elusters. Numerieal simulations suggest that shocks with the Maeh number of M = 3-5 
arise in late phases (i.e. after eore erossing) of cluster mergers (e.g. Ricker and Sarazin 2001). Sueh 
shoeks, eharacterized by significant density/temperature jumps and sometimes with radio emissions 
(e.g. Feretti 2012), have been studied intensively with X-Ray and radio observations. It is, however, 
not easy to determine the eollision parameters in late phase mergers, sinee original ICM temperature 
and density before shoeks passage is not elear. An example is the nearby late phase merging clusters 
Abell 3667; it exhibits clear evidence for shock heating of the ICM (e.g. Briel et al. 2004, Nakazawa 
et al. 2009, Finoguenov et al. 2010, Akamatsu et al. 2013) together with two bright diffuse radio 
emission, namely relics (e.g. Roettiger et al 1999), but its merger geometry is still unelear. In ease of 
another famous late phase merger, the Bullet eluster (1E0657-56), the merger geometry looks elear 
(e.g. Markeviteh et al. 2002) although the original ICM temperature itself is not well determined. 

We can more easily understand the merger geometry in early phase (i.e before core erossing) 
merging clusters, such as Abell 222-223 (e.g. Werner et al. 2008), Abell 399-401 (e.g. Fujita et 
al. 1996, Fujita et al. 2008) pairs although they do not show strong evidence for shoek heating yet. 
Cygnus-A eluster of galaxies, for example, is one of a few pre-core erossing nearby merger with 
evidenee for shook heating in between the two oores (Markeviteh et al. 1999, Sarazin et al. 2013). 
In this ease, there is no clear deteotion of surfaoe brightness jump, and the shook looation is rather 
ambigious. 

In the present work, we report on Suzaku and XMM-Newton observations of an early phase 
merger CIZA J1358.9-4750 whieh elearly exhibits shook signature. We use the Hubble constant Hq 
= 67 km s~^ Mpc“^, together with the cosmological parameters Qm = 0.31 and (Ia = 0.69 (Planok 
Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014). The errors refer to 90% oonfidenoe for a 
single parameter. 

2. CIZA J1358.9-4750 

This object, CIZA J1358.9-4750, is listed in the CIZA oatalog (Ebeling et al. 2002, Kocevski 
et al. 2007), whieh colleots X-Ray selected elusters of galaxies in the Zone of Avoidanoe (typioally 
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Fig. 1. A 0.6-10 keV X-Ray image of CIZA J1 358.9-4750 taken with Suzaku XIS3, shown without background 
subtraction. Spectra are extracted from regions indicated in green. Light-blue circles indicate the 3 point sources 
discussed in section 3.1. 

|b| < 20, where b is the Galactic latitude). By inspecting ROSAT and XMM-Newton short exposure 
images, 3.2 ks and 4.8 ks respectively, we identified this object as a nearby (redshift of ~ 0.07, that is, 
angular diameter distance of 290 Mpc) major merger candidate. The 0.1-2.4 keV ROSAT flux reaches 
1.92 X 10“^^ ergs cm“^ s“^ which is the sixth highest in the CIZA catalog excluding point source 
candidates. It has two distinct X-Ray humps in the north-west (NW) and south-east (SE) directions, 
separated by ~ 14 arcmin (approximately 1.2 Mpc), which are both considered to represent clusters of 
galaxies. Furthermore, a clear X-Ray enhancement is seen to connect them. A bright elliptical galaxy, 
considered to be the cD galaxy, is located right at the center (within 1 arcmin = 84 kpc) of each X- 
Ray hump. The two galaxies (2MASX J13590381-4751311 and 2MASX J13581085-4741243) at the 
center of the SE and NW clusters have optical redshifts of 0.0745 and 0.0709, respectively (from NED 
database). 

3. SUZAKU Observation 

With Suzaku, we observed CIZA J1358.9-4750 on 2013 January 21-23, covering the two X- 
Ray humps. The XIS was operated in the normal full-frame clocking mode. In this paper, we employ 
the data of both 3x3 and 5x5 edit modes acquired with XISO and XIS3. With the standard data 
screening, we obtained an exposure of 61.7 ksec for the XIS. The 0.6-10 keV image of the XIS3 
detector is shown in figure 1. It clearly reconfirms the two X-Ray humps and the bridge-like X-Ray 
emission between them. 
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3.1. Spectral analysis 

As shown in figure 1, we defined 24 regions shown in green for our speetral analysis. Among 
them, eireles labeled N and S are for the NW and SE elusters, respeetively, with a radius 2M7 
(180 kpe) eentered on the X-Ray peaks. The regions B1-B6 and C1-C6 eover the eentral parts (near 
the suggested eollision axis) of the bridge, augmented by the series A1-A6 and D1-D4. Eaeh region 
has a size of 4'.0 x 1'.5 or 3'.0 (335 kpe x 130 or 260 kpe). Then, using XISO and XIS3, on-souree 
speetra were aeeumulated over the 24 regions individually. The XISO and XIS3 data were added to¬ 
gether, and instrumental baekground (non X-Ray baekground: NXB) was estimated from night Earth 
data by xisnxbgen (Tawa et al. 2008). The baekground was then subtraeted from the speetra of eaeh 
region. Note that beeause the regions A5 and A6 happen to be on the XISO segment-A, whieh is not 
well working, we only used XIS3 data. 

We eonsider three sky baekground eomponents; the loeal hot bubble (EHB) emission with a 
temperature of 0.1 keV, that from the Milkyway halo (MWH) with 0.3 keV, and the eosmie X-Ray 
baekground (CXB). The temperatures of EHB and MHW were determined from an off-souree Suzaku 
speetrum, obtained in an observation of the star HD 125599, whieh is loeated 4° away from CIZA 
J1358.9-4750 at a similar Galaetie latitude. The normalizations of EHB and MHW were determined 
by fitting all the 24 extraeted speetra of CIZA J1358.9-4750 simultaneously, assuming EHB and 
MWH have a eonstant surfaee brightnesses in the observed regions. The 2-10 keV surfaee brightness 
of the CXB was assumed as 6.38 x 10“® erg em~^ sr“^ using the parameters of Kushino et al. 
2002, and its photon index was fixed at T = 1.41. All these parameters are listed in table 1. 

Eor speetral fits response matrix files and aneillrary response files were ealeulated using xis- 
rmfgen and xissimarfgen, respeetively (Ishisaki et al. 2007). By fixing the parameters of the baek¬ 
ground eomponents as determined above, the 24 NXB-subtraeted speetra in energy band of 0.6- 
10 keV were fitted individually by a model of apeo[EHB] -i- wabs x (apeo[MHW] -i- powerlaw[CXB] 
-I- apeo[ICM]) [model in XSPEC]. The interstellar absorption was represented by the wabs faetor, 
where in hydrogen eolumn density fixed at Nh = 11.5 x 10^° em“^ as taken from the 21em mea- 
surments (Kalberla, R M. W. et al. 2005). We left free the temperature, abundanee, redshift, and 
norm parameters of the ICM, but the abundanee and redshift parameters were fixed as 0.2 and 0.07, 
respeetively, in low statistie speetra, i.e, having eount rates less than 2.0 x 10“^ ets s“^. An energy 
band of 1.7-1.9 keV was ignored beeause of the prominent instrumental feature. The speetra of S, N, 
A3 and A4, B3 and B4 are shown in Eig.2. The 24 fits were all sueeessful, with x^/d.o.f in the range 
of 0.79 to 1.18, and yielded the best-fit parameters as given in table 2. Thus, the ICM temperatures of 
the individual regions have been determined with a typieal 90% error of ±1 keV. 

There are 3 point sourees with the 0.4-7.2 keV X-Ray flux higher than 1.0 x 
10“^"^ ergs em“^ s“^, estimated from the XMM-Newton short exposure data, as detailed in seetion 
4. By exeluding regions 1'.5 around these sourees, we estimated their eontamination effeets to our 
Suzaku temperature analysis. Only the result from region D2 ehanged signifieantly (~ 1 keV), while 
the effeet in other regions is negligible, i.e. less than 0.1 keV. Beeause we have eut the 3 point sourees. 


Table 1. Summary of the fixed prameters for background components. 



LHB 

MHW 

CXB 

kT{ keV)/r 

0.10 

0.29 

(1.41) 

norm (x 10“^)’!' 

35.61 

2.59 

(0.96) 


^ : nomalization of apec and powerlaw models in XSPEC verl2.8.0 scaled with factor 1 /4007r. 



Fig. 2. Suzaku (XIS0+XIS3) spectra of CIZA J1358.9-4750, extracted from regions labeled S, N, A3, A4, B3, and 
B4 in figure. Red, cyan, dark blue and magenta indicate the ICM, LHB, MHW, and CXB components, respectively. 


the CXB normalization could be smaller than those provided by Kushino et al. 2002. Thus, we also 
checked if the results can change with the smaller CXB normalization derived from Lockman hole 
observation with almost no source contamination (5.7 x 10“® erg cm~^ sr“^, Suzaku observation 
ID 101002010; see Nakazawa et al. 2009) and found it to be insignificant, i.e. less than 0.1 keV. From 
the above, we replaced the results for region D2 using the source-excluded spectra. We also checked 
the effect of fixing Nh through-out this analysis, and also found that the difference cause by making 
the value free is negligible (less than ~ 0.4 keV) to the spectral parameters of the ICM component. 

3.2. Temperature distribtuion 

The derived temperature distribution is plotted in figure 3. Center of the SE and NW clusters 
were found to have a temperature of 5.6 ± 0.2 keV and 4.6 ± 0.2 keV, respectively. The small differ¬ 
ence, approximately 20%, indicates that they have similar mass from M -T relation (Vikhlinin et al. 
2006). These temperature suggest that Virial radius, i.e. r 2 oo, of the SE and NW clusters are 1.7 Mpc 
and 1.5 Mpc, respectively. Using /9-model assuming (3 = 0.6 and = 180 kpc at 0.5-10 keV, lumi- 
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Table 2. The parameters of Temperature, Abundance, normalization. The letters specify the regions shown in hgure 1. 



At 

A2 

A3 

A4 

A5 

A6 

kT( keV) 

6.5o;°:?3 

‘ -04-0.84 

9-26;f°I 

O.OO_0 88 

0 09+1-97 

o.^^_l 30 

K 7(4+1-06 

0. f U-0.88 

Z0 

0.38t°;i« 

O.ZZ-o 11 

fixed 

fixed 

n 90+0.26 

U.25-0.20 

fixed 

X^/d.o.f 

137/147 

212/205 

130/112 

88/100 

108/115 

69/77 


B1 

B2 

B3 

B4 

B5 

B6 

kT( keV) 

7 +0-74 

^ ■'^-‘^-0.59 

8-36t°;®^ 

9.14;/" 

a 07+0-75 
-0.61 

7 1 7+0.79 
' -0.70 

ft 07+0.58 
_o.54 

Z0 

0 Tfi+0.12 

0.46t°;}t 

0.25;/" 

o'-^y-o.io 

n 10+0-12 
o’- 4 y_o.io 

n 1 Q+0.10 

X^/d.o.f 

267/248 

235/249 

244/233 

222/210 

218/211 

183/217 


Cl 

C2 

C3 

C4 

C5 

C6 

kT( keV) 



8.20;;" 

Q-QQ-ofi 

6.6o;°:°° 

6 11 +109 
O'-Li-o.go 

Z0 

0.30l°;i° 


0.27;/" 

0.18;°;" 

0.25;/" 

0 11+0.12 

X^/d.o.f 

262/245 

207/224 

192/205 

178/186 

160/172 

167/165 


D1 

D2 

D3 

D4 

s 

N 

kT( keV) 

0.00-0.87 

5.06tJ:°^ 

6.26;;;® 

5.89;°:?° 

p; fin+o io 

A p:c:+0.16 
^■^0-0.16 

Z0 

fixed 

fixed 

n 10+016 
O-4y-0.09 

0 21 +0'i6 

Q QO+O-Oe 

o.oy-0.02 

'9-^^-0.04 

X^/d.o.f 

89/112 

70/72 

159/160 

126/131 

795/800 

633/629 


nosity within r 2 oo of the former cluster is estimated to be 2.6 x 10^^ erg s“^, scaled from the central 
180 kpc flux. Similarly, those of the NW cluster is estimated to be 1.8 x erg s“^. To search for a 
signature of a cooling core, we extracted spectrum of N and S within a radius of 90 kpc (1'.08), and 
derived temperatures of 5.7 ± 0.3 and 4.4 ± 0.3, respectively. Therefore, neither of them show strong 
evidence for a cool core. 

Compared to the two clusters, A3, B3 regions in the bridge region show higher temperatures 
typically by a factor of ~1.5. The hot region extends to almost perpendicular to the axis connecting 
the two cores, namely, A3, B3, and C3 regions. The westward regions D1-D4 appear relatively 
isothermal. The measured temperatures exhibit a particularly significant increase from ~ 7 keV in the 
region A4-B4-C4 to ~ 9 keV in A3-B3-C3. 

3.3. Entropy distribution 

Sharp temperature increase is suggestive of the presence of a shock discontinuity. To identify 
whether this is caused by shock or adiabatic compression, entropy distribution are plotted in figure 4. 
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Fig. 3. Temperature distributions of CIZA J1358.9-4750 in the A1-A6 and D1-D6 regions (left), and B1-B6 and 
C1-C6 (right). Abscissa represents the angular distance along the bridge axis. Note that statistics of the A5 and A6 
regions are about a half, since we could not add data from XISO. 

For simplicity, we employ the “astrophysical entropy” expressed as K = kT/n^^^. Since it preserves 
heat input to the ICM, entropy is a good indieator of shocks. Here, we simply estimated the gas 
density by assuming a eylindrical geometry with a radius of 335 kpe for regions B and C which is the 
width of these regions. In regions A and D, more simpler cuboid geometry with line-of-sight depth 
of 670 kpc and horizontal width of 335 kpe x 130 or 250 kpe is adopted, because the over-all shape 
of the eluster around the periphery is not resolved. At regions of A4-B4-C4 to A3-B3-C3, the same 
place where temperature increase rapidly, the entropy shows a elear jump of a factor of ~ 1.4. Thus, 
temperature increase is not only produced by adiabatic compression but also shock. There are few 
example of early merging clusters with clear evidence for shock. In addition, we also note that the 
shoek is detected in relatively bright region and is easy to observe in detail thanks to its brightness. 
Therefore, this objeet is important to understand evolution of shoeks as well as nature of shook. 


4. Flux jump in the XMM-Newton image 


Shoek is found from Suzaku entropy distribution, but further information about its detailed 
orientation is hampered by the limited angular resolution. To seareh for a possible surface bright¬ 
ness jump eaused by the shoek suggested by figure 4, we further utilized an archival XMM-Newton 
snapshot data (ObsID = 0204710501; aoquired in 2004 February 21). Basie reduction and calibration 
of the European Photon Imaging Camera (EPIC) data were earried out with SAS (ver. 13.5). In the 
screening proeess we set FLAG = 0, and kept events with PATTERNS 0-12 for the MOS eameras and 
those with PATTERNS 0-4 for the PN camera. By examining lightcurves extraeted in 10.0 — 14.0 keV 
and 1.0 — 5.0 keV from souree free regions, we disearded time intervals affeeted by hard-band and 
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Fig. 4. Entropy distributions of CIZAJl358.9-4750 in the A1-A6 and B1-B6 regions (left), and C1-C6 and D1-D6 
(right). Abscissa represents the angular distance along the bridge axis. 

soft-band flares, respectively. The obtained MOS and PN exposures are 4.1 ks and 1.8 ks, respectively. 

The raw event image in the 0.5-4.5 keV band is shown in Fig.5. We excluded 4 point sources 
selected with the same criteria as we did for the Suzaku spectral analysis. We then rotated the image 
by 25 degree counter-clockwise and produced brightness profiles to find a clear jump at the position 
1500 arcsec in the figure. East region has biggest intensity jump and west is weaker. Although the 
XMM-Newton exposure is too short to determine the ICM temperature, the location of the brightness 
jump is right within the Suzaku high-temperature regions. 

The shock front detected in the XMM-image is a bit inclined from the axis between the two 
peaks. While the bridge axis is 43 degree clockwise to the north, the shock normal is inclined by 
25 degree as noted above. Because the region defined in Fig.5(a) is aligned to the former, we also 
checked the temperature profile with a slightly rotated regions following the later angle, to find virtu¬ 
ally the same results; from SE to NW, it showed a temperature jump from S-filgira keV to T.hloJ® keV 
in regions corresponding to “mid”. Note that both regions have a size of 1'.5 x 3'.0. 

5. Discussion 

5.1. Merger phase 

Using Suzaku, we confirmed that the X-Ray image of CIZA J1358.9-4750 has two humps, 
corresponding to the two clusters without cool cores, and a bridge structure. In addition, we found 
that there is a temperature jump within the bridge. This X-Ray image and temperature distribution are 
similar to those seen in numerical simulations of an early phase of cluster major merger (e.g Akahori 
& Yoshikawa 2010; Takizawa 2008). The good positional coincidence between the cD galaxies and 
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Fig. 5. (a) A 0.5-4.5 keV image of CIZA J1358.9-4750 with XMM-Newton. Surface brightness profiles are ex¬ 
tracted from regions indicated in white, red and green, respectively. Small light-blue and yellow circles are the 
excluded regions around point sources with radius of 20" and 30", respectively, (b) 0.5-4.5 keV surface brightness 
distributions of CIZA J1358.9-4750. Horizontal axis is in arcsec with arbitral zero point. 

the X-Ray peaks provides another support to this view, where in the late stage of merger the coinci- 
denee tends to be not observed (e.g., Abell 2163; Okabe et al. 2011). 

As shown in Fig.5(b), the X-Ray surface brightness profiles derived with XMM-Newton ex¬ 
hibit a clear jump at the position where the ICM temperature measured with Suzaku exhibits a steep 
change. This good positional coincidence implies a pressure jump in that region, that is, the presence 
of a shock discontinuity. Thus we can apply the Rankine-Hugoniot condition to this jump. Employing 
the largest temperature increase, 6.97lo;gf keV in B4 to 9.14lJ;22 keV in B3, we obtain a Mach num¬ 
ber of 1.32 ± 0.22. Since the sound velocity in the ICM is ~ 1360 kms“^ at 6.97 keV, i.e, at the 
pre-shock region, the colliding velocity is approximately ~1800 kms~^. Because the width of the 
spectral fitting regions are as narrow as 1'.5 which is smaller than the Suzaku PSF of 2' half power 
diameter, the spectral variation shall be smeared. The true temperature jump therefore could be a bit 
higher. In other words the shock Mach number presented here is lower limit. 

In the present study we fitted the ICM emission with only one APEC model (IT), because 
it provided acceptable fits. According to numerical simulations (e.g., Takizawa 2008; Akahori, 
Yoshikawa 2010), it is plausible that there are multi-temperature structures in the projection along 
the line of sight toward the shock-heated regions (A3-B3-C3), i.e. the shock-heated very hot com¬ 
ponent and the ambient, original component. Therefore, our assumption of single temperature may 
provide a lower limit on the temperature increase, and the Mach number and the colliding velocity 
could be higher than as quoted above. 

Numerical simulation of merging clusters suggest that there should exist a forward-shock 
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and reverse-shock as shown schematically in fig.5 of the review paper presented by Markevitch & 
Vikhlinin (2007). In this observation, clear symptom of reverse-shock was not discovered. This 
could be merely by the lack of angular resolution (Suzaku), exposure (XMM-Newton) and strength 
of reverse-shock. But, in Fig.5(b), east region exhibit brightness jump slightly at 120" = 2' (170 kpc) 
from the identified discontinuity, resembling the reverse shock. With Mach 1.32 and the sound ve¬ 
locity of 1360 km s“^ in the ICM, post-shock velocity can be derived to be 1200 km s“^. If we 
assume a 1:1 merger and the post shock region has no bulk motion as the simplest approximation, 
the shock-front velocity on the sky plane becomes 1200 km s“^ and the reverse shock has the same 
velocity from symmetry. Then the narrow width of 170 kpc suggests that the shock is young, i.e. only 
~700 Myr has passed since its birth. 

5.2. Merger geometry and its implication 

Because the shock front is clearly visible, the collision axis is considered to be nearly parallel 
to the sky plane. To confirm this, we utilized the NED database to identify the optical redshifts of the 
two clusters. We found 8 and 11 redshift-measured galaxies within 455 kpc of the SE and NW clusters, 
respectively. The average redshifts are 0.0721 and 0.0737 in the NW and SW humps, respectively. 
Since the difference, A v = 480 kms“^, is significantly smaller than the estimated colliding velocity, 
the two clusters are confirmed to be almost on the sky plane. 

One open question is the fact that the (candidate) cD galaxies have significantly different 
redshifts, with At = 1000 km s“^. However, we do not have any additional information on this issue 
and leave it for future investigations. 

6. Summary and Conclusion 

We found that CIZA J 1358.9-4750 is made of two medium sized clusters separated by 1.2 Mpc 
which are connected by a bridge of enhanced X-Ray emission. The bridge region is hotter than the 
main emission of each cluster and the position of temperature increase coincides with a brightness 
jump. We identified this structure with a shock fronts, and derived a Mach number of 1.32 ± 0.22. 
Thus, CIZA J1358.9-4750 is concluded to be a precious example of early-phase major merger with 
clear evidence of shock and shock heating, produced as the two clusters of similar sizes are making a 
head on collision along the sky plane. 
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